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Abstract. In a couple of articles (Ford G W and O'Connell R F 1996 Nature 380 
113 and 2002 J. Phys. A: Math. Gen. 35 4183) it was argued that the standard result 
for the derivative of the hyperbolic cotangent in the literature, dcothy/dy = — csch 2 ?/ 
is incomplete and the correct expression should have an additional term proportional 
to the Dirac delta function. The purpose of this paper is to demonstrate that this 
claim is incorrect. 
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In [TJ Ford and O'Connell published the formula 

— cothy = — csch 2 ?/ + 2 5(y) (1) 

for j/El. Their argument showing that it is correct can be summarized as follows. A 
function f(y) = cothy, satisfying f(—y) = —f(y), is identically rewritten as 

f(y) = sign(y) {l + g(\y\)} (2) 

where an auxiliary function g(y) = 2/ (e 2y — 1), and sign(y) = ±1 depending on the sign 
of y. It was then claimed that while the derivative of the second term, sign(j/) g(\y\) 
equals — csch 2 ?/ for every y 0, the derivative of the first term, sign(y) • 1 yields 28(y). 
Consequently, the derivative of f(y) should contain the extra delta-function term. It 
is clear that this logic is faulty. Indeed, we could repeat the same argument for an 
arbitrary odd function f(y) with a constant c ^ in place of factor 1 of the first term, 
which leads to the conclusion that df/dy would contain a term 2cS(y). The error is, of 
course, in disregarding the delta-function term, — 25 (y) stemming from the derivative 
of the second term (for its detail, refer to equation (48) and the discussion thereafter in 
[2]). In fact, by means of a partial fraction series for the hyperbolic cotangent [3] 

00 . 

cothy = V — — — , (3) 

y + iku 

k=—oo 

we can explicitly show that dcothy/dy = — csch 2 ?/. Here, we used [3] 

1 °° 1 

cs ° 2y = shr^ = (y - k71 ) 2 (4) 

with sinhy = —ismiy. 

The derivation of equation (pp) was motivated by the need to verify consistency 
between a quantum statistical mechanical quantity and its classical counterpart to 
be identical to the expression obtained from the quantum-mechanical quantity in the 
limit of h — > 0. However, equation (3) in [1] (and also equation (2.11) in [4J) for the 
autocorrelation of the quantum-mechanical random force F(t) in the case of constant 
friction £, 

- (P(t) F(f) + F(t') F(t)) = - [ duohu coth(Huj/2kT) cos[uj(t - t')} (5) 

2 \ / 7T J 

= kT( — coth[7r kT(t - t')/H] (6) 

(JjL 

is actually incomplete while equation (4) in pQ (and also equation (2.12) in [4]) for its 
classical counterpart is correctly given as (F(t)F(t')) = 2kT(5(t — f). The integrand in 
equation j5]) is an even function with respect to the integration variable oo. Hence, the 
integral (understood in the distributional sense) may be performed from u = — oo to 
+oo and can be evaluated in closed form with the aid of identity ([3|). We then obtain 
an additional term to j6]) for the autocorrelation, arising exactly from k = in (j3J), 
2kT( 



71 



poo 

/ du cos u(t-t') = 2kT(5(t — t') . (7) 
Jo 
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This precisely corresponds to the classical autocorrelation while equation ((6)) represents 
quantum fluctuations for the autocorrelation about the classical mean value j7]) so that 
(Ej) vanishes in the limit of h — > 0. It was, however, argued in [2] (c/. equations (44) 
and (45) therein) that simply from equations (DQ) and — csch 2 y can be interpreted as 
the portion of dcoihy/dy representing quantum fluctuations about the classical mean 
value 25 (y) such that dy csch 2 y = 0. This claim is obviously incorrect since the 
additional term (13) was just omitted (note also that csch 2 y is an even function anyway). 
The expression rfcoth y/dy has no direct physical meaning whatsoever. 

In [5] Ford and O'Connell gave some additional detail on the "derivation" of formula 
(PQ). Motivated by the observation that cothy increases by +2 as y goes from — oo to 
+oo, they identified cothy to Re{coth(y + ie)} with e — > where 

. , , sinh2?/ . sin2e 

coth(y + ie) = — i — , (8) 

cosh 2y — cos 2e cosh 2y — cos 2e 

and argued that Re {coth(|/ + ie)} can be understood as a distribution to yield the 
additional term 25(y) for <icoth /dy. However, it is incorrect. To show this, we consider 
the identity obtained from the interplay between the distributions and the theory of 
moments [6] 

n=0 

where //„ = J™ dy g(y) y n are the moments of g(y). For g(y) = Im {coth(y + ie)}, we 
have H2m+i = and 

2 (_ 1 )m+l 7r 2m+l 

^2m — ~ — -D2m+1 (10) 

2m + 1 

where B n are the Bernoulli numbers. Here we used [3] 

x 2m 2(2m)\ ^sin2kaTT , . 

dx — : = ~ ~ -> 11 

cosh x — cos 2a,7T sm 2an z — ' ^m+i 

k=l 



for o ^ 1/2, and 



oo 



sin2£;a7r (-l) n+1 (27r) 2n+1 



^ fc 2n+1 2(2n + l) 

fc=i v ' ; 



B 2n+ i{a) (12) 



where B n (a) are the Bernoulli polynomials with the property, B n (0) = B n . From 
and (El) we then obtain 



coth(y + ie) = Re {coth(y + ze)} + z 1 (2 ^ , (2m + ^ ^ (2m) (l/) • (13) 

The first term on the right hand side corresponds to the principal part V {coth(?/ + ie)} 
only, which vanishes at y = 0, while the second term consists of the distributional 
contributions as for the well-known identity l/(y ± ie) = V(l/y) =F inS(y) where 
V(l/y) = y/(y 2 + e 2 ). As a result, cothy cannot be identified with Re{coth(y + ie)}, 
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nor can Re{coth(y + ie)} be understood as a distribution. Further, by using identity 
([3j) with y + ie in place of y <G R we can get 



and thus cothy — V {coth(y + ie)} = 1/y — Vil/y). This demonstrates that coth?/ and 
its derivative cannot be distributions either as none of 1/y, V(l/y),V {coth(y + ie)} and 
their derivatives can. In addition, we see from identity ([3D that coth(±oo) = ±1^0, 
as compared to 1/y = at ?/ = ±oo, stem simply from a sum of all terms except 
for k = 0. Therefore, the increase of cothy by +2 as y goes from — oo to oo is not 
ascribed to the additional delta-function term in (Op) which should be removed. Actually, 
L(y) := cothy — 1/y is called the Langevin function and plays an important role in 
quantum statistical mechanics [7]. 

In [2] how singular functions can define distributions was discussed by regularizing 
the singularities in the standard way, namely, by either analytical continuation or 
extraction of the Hadamard finite part. Special attention was paid to the distribution 
(csch 2 y) st obtained from the standard definitions. It was then shown that d cothy / dy = 
— (csch 2 y) st with no additional delta-function term. It was also argued that the classical 
mean value ([7]) can be extracted from equation ([0]) or Xg(X(t—t')) = —kTQ X (csch 2 [A(t— 
t')]) st with A = nkT/h in the limit of A — > oo. Using the moment asymptotic expansion 
([9j), Xg(X(t — £'))|a-+oo exactly reduces to the classical value —kT(fi 5(t — t') with 
A*o = X^o( csc h 2 !/ — l/?/ 2 ) dy = —2 (refer to equations (49)- (53) in [2]). However, it may 
be argued, from the physical viewpoint, that equation ([6j) at t = t' cannot reflect the 
semiclassical limit h — > 0, corresponding to t — tf ^> 1, in which identically vanishes. 
As a result, this distributional concept is not needed to sustain the consistency between 
classical and quantum autocorrelation in the limit of h — > when the additional term 
(I7|) is correctly considered. 

The author is grateful to S. Winitzki (Munich) for stimulating discussions. 
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